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Although the sex ratio of Odonata at emergence has received much attention, we are still far from under- 
standing the exact causes of its variability and imbalance. In this paper we studied the sex ratios at emergence 
in natural populations of two Gomphus species based on samples of exuviae taken from two European 
lowland rivers. We hypothesized a possible relationship between the water temperature during larval devel- 
opment and the sex ratio at emergence. Sex ratio records exhibited no consistent bias towards one sex but 
varied between habitats and years in both species. We found correlations between sex ratio and water 
temperature in the year preceding emergence. Furthermore, the correlation between sex ratios and water 
temperature was in opposite directions in the two congeneric species, which may be attributed to dif- 
ferences in their voltinisms. We conclude that the effect of water temperature can be mediated through 
cohort-splitting; temperature-dependent development of minor cohorts, including unequal proportions of 
males and females due to the faster development of male larvae, affects the sex ratio at emergence. The 
supposed effect does not cause a long-term consistent bias, but may explain the year-to-year variations. 


Keywords: dragonfly; Gomphus flavipes; Gomphus vulgatissimus; Rivers Tisza and Szamos; voltinism; 
cohort-splitting; exuviae 


Introduction 


The sex ratio of Odonata at emergence has received much attention. Several reviews, summarizing 
all records available (Corbet, 1983, 1999; Corbet & Hoess, 1998; Lawton, 1972), are consistent 
in the main conclusions that sex ratios at emergence often exhibit an imbalance; and that this 
imbalance indicates a slight excess of females in the Anisoptera and the opposite in the Zygoptera. 
However, besides this general tendency, sex ratios can differ greatly among taxa, even within a 
species (Corbet & Hoess, 1998). 

It has long been known that sex determination in Odonata in almost all cases is of the XO/XX 
type, the male being the heterogametic sex (Kiauta, 1969). This mechanism produces a sex ratio 
of unity in the zygote. Any variability and imbalance in sex ratios observed at emergence is 
therefore related to the egg and/or larval stage (Corbet & Hoess, 1998). Several assumptions and 
hypotheses have so far been described to explain the biased sex ratios (Baker, Forbes, & Proctor, 
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1992; Corbet, 1957; Fuselier, Decker, Lunski, Mastel, & Skolness, 2007; Gribbin & Thompson, 
1991; Johansson, Crowley, & Brodin, 2005; Lawton, 1972; Purse & Thompson, 2003; Wong- 
Muñoz, Cérdoba-Aguilar, Cueva del Castillo, Serrano-Meneses, & Payne, 2011). Yet, as far as 
we know, there is little evidence supporting any of these hypotheses. Hence, the causes of sex-ratio 
imbalance has remained largely unknown, as has its biological significance (see Corbet, 1999). 

The majority of previous studies on sex ratios at emergence focused on a simple scenario 
in which the sex ratio exhibits a consistent bias towards one sex in a given species (Crowley 
& Johansson, 2002; Falck & Johansson, 2000; Fuselier et al., 2007; Johansson et al., 2005; 
Wong-Mufioz et al., 2011). However, many species show widely varying sex ratios between sites 
and years without a consistent bias (e.g. Corbet & Hoess, 1998; Miiller, 1995). This suggests 
that environmental factors may be influential in altering the sex ratios, although the correlation 
between them is yet to be identified (Corbet & Hoess, 1998). 

In temperate regions, besides photoperiod, water temperature also plays a major role in larval 
development, affecting larval growth rate, seasonal regulation and voltinism (Corbet, 1999). At 
a given latitude, exposed to the same photoperiod, differences in larval development are mainly 
attributed to different temperature regimes (cf. Miiller et al., 2000). Since the variability in the 
sex ratio is most likely to be related to the larval stage (see above), we hypothesize a possi- 
ble relationship between the water temperature during larval development and the sex ratio at 
emergence. 

To test this hypothesis in this paper we studied natural populations of two riverine dragonfly 
species (Odonata: Gomphidae), Gomphus flavipes flavipes (Charpentier, 1825) and Gomphus 
vulgatissimus (Linnaeus, 1758). To our knowledge, it was the first attempt to investigate the 
relationship between sex ratio and water temperature. 


Materials and methods 


Sampling sites and data collection 


Sex ratio records of G. flavipes and G. vulgatissimus were based on the systematic collections 
of exuviae. Collections were carried out at four different sites along the Tisza river system in 
Northeast Hungary: (1) along the upper reach of the River Tisza at Jand, upstream of the mouth 
of River Szamos; (2) along the upper reach of the River Tisza at Vásárosnamény, downstream of 
the mouth of River Szamos; (3) along the River Szamos at Olcsva, near the confluence; and (4) 
along the dammed middle reach of the River Tisza between Tiszacsege and Tiszafiired (henceforth 
Tiszafüred) (Table 1). The three sites close to each other (1-3) were sampled in the same year 
(from late April to mid August in 2008), while the dammed site (4), much further downstream, was 
sampled in four different years (from late April/early May to mid August in the years 2001-2003 
and 2009). The study sites represented a range of hydrological and hydromorphological conditions. 
Furthermore, both different sites and years were exposed to various temperature regimes. 

The collections of exuviae were made on standardized areas, 20 m long stretches of the river 
banks. The stretches at a given sampling site differed in river bank characteristics ranging from 
almost flat to a steep bank normally covered with sandy sediment, but with paving at certain 
stretches. Riparian vegetation also varied from dense to scattered stands of trees or bushes (mainly 
Populus and Salix) with different undergrowth cover (especially Rubus). During collections the 
riverside ground and the vegetation were searched extensively for the exuviae in an area of the 
bank 6 m in width. All the exuviae found were collected and later counted, identified and sexed 
for each day. 

Samples were taken regularly and frequently during the whole emergence period (see Farkas, 
Jakab, Tóth, Kalmar, & Dévai, 2012); daily for the most part of the emergence period at sites 1-3, 


Table 1. Characterization of the sampling sites and collection method. 











River Discharge Tinn: Tyan Tyuly Collection Sampling Collection 
Collection site Geo-coordinates River km width (m) (m3s—!) (°C) (°C) (°C) season unit (m) frequency 
Tiszafüred, R. Tisza, 1.b. 47°40'N 20°50'E 433—450 ~175 574 12 1.2 23.3 2001-2003, 2009 6 x 20 Twice per week 
Jánd, R. Tisza, r.b. 48°06'N 22°21'E 687-689 ~117 249 10 0.8 20.3 2008 3 x 20 Daily* 
Vásárosnamény, R. Tisza, 1.b. 48°07/N 22°20/E 684-685 ~121 392 10.5 0.9 21.3 2008 3 x 20 Daily* 
Olcsva, R. Szamos, 1.b. 48°05'N 22°20'E 3—4 ~62 147 11 0.9 22 2008 3 x 20 Daily* 
Notes: 1.b.: left bank; r.b.: right bank; Tann.: mean annual water temperature; Tyan: mean water temperature in January; Tyuy: mean water temperature in July (based on 20-year data). 
*When emergence intensity decreased and seasonal floods moderated, the stretches were checked twice a week (in the last two weeks of the emergence at most). 
Table 2. Sex ratios of the studied species at emergence. 
Collection LS LS VS LS 
Species Collection site season Onset n LS EMo-25 EMo-35 EM75_100 EMoo_100 
G. flavipes Tiszafüred 2001 27 May 1888 49.2 56.23 54.52 44.61 46.15 
2002 25 May 3105 48.0 54.37 51.70 44.49 37.76 
2003 27 May 4883 51.2 57.39 57.38 44.92 44.71 
2009 25 May 3145 48.3 57.02 54.57 41.76 39.33 
Jand 2008 10 June 54 51.9 57.14 51.36 46.15 64.29 
Vásárosnamény 2008 1 June 156 53:2 63.02 61.89 52.83 59.75 
Olcsva 2008 31 May 183 50.8 57.14 57.18 50.00 58.44 
G. vulgatissimus Tiszafüred 2001 4 May 250 48.8 62.71 62.83 36.07 28.25 
2002 4 May 298 52.0 62.55 60.80 43.57 42.40 
2003 10 May 209 44.0 48.77 47.75 38.16 33.30 
2009 28 April 120 50.8 59.72 58.44 39.69 40.71 
Jand 2008 12 May 66 43.9 78.24 62.75 12.50 13.79 
Vásárosnamény 2008 6 May 88 43.2 40.40 42.34 37.50 61.36 
Olcsva 2008 3 May 91 38.5 65.16 54.36 16.75 11.52 


Notes: Onset: onset of emergence; n: total number of exuviae that can be sexed; %g': percentage of males of the total; %8 EMo-25: percentage of males of the total from the beginning of emergence until EM25 was 
reached; %c" EMo-35: percentage of males of the total from the beginning of emergence until EM35 was reached; %c" EM75_100: percentage of males of the total between EM75 and EMj09; %8 EMop-_100: percentage 


of males of the total between EMog and EM jo0. 
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where the water level often changes, and twice a week at the dammed site, where the water level 
is almost constant due to the barrage. 

Sampling methods met the criteria listed by Corbet and Hoess (1998) for ensuring that sex ratio 
data are accurate and comparable. Due to the low densities, though, the sample size sometimes 
did not exceed the recommended minimum of 99 exuviae. We considered this in our analyses (see 
below). 

Sex ratio was expressed as a percentage of males (%<") of the total number of exuviae collected 
during the whole emergence period. Percentages of males were also calculated for distinct periods 
of emergence (EM): %c" from the beginning of emergence until EM25 and EM35 was reached 
(%8 EMo-25 and WF 'EMo-35 respectively), as well as %g between EM7s5 and EM jo and between 
EMoo and EM j99 (%g EM75_1 09 and %o" EMop_109 respectively). The former represent distribution 
of the sexes during early phases of emergence, while the latter represent the late phases. 

Long-term (1989-2009) water temperature data for the sampling sites were obtained from the 
Environmental and Water Management Directorate of the Tiszántúl Region and from an internet 
database (http://www. vizadat.hu). Water temperature data were derived from the nearest measur- 
ing stations and were measured daily at 7:00 CET by the local Water Management Directorates. 

In order to study the effect of water temperature on sex ratios, the mean temperatures of two 
distinct periods of larval development were used. At our sampling sites (at 47-48°N in two 
medium-sized rivers), life cycles of 2 to 3 years long are suggested for both species (Braune, 
Richter, S6ndgerath, & Suhling, 2008; Müller, 1995; Müller et al., 2000; Sondgerath, Rummland, 
& Suhling, 2012). Therefore, water temperatures from two years preceding emergence were taken 
into account. This time window should cover the major part of the larval stage. However, the period 
when larval growth is interrupted (from diapause induction until spring when rising temperatures 
allow further development) should be excluded. Due to the fact that day length and temperature 
thresholds for diapause are variable in Odonata (Corbet, 1999; Flenner, Richter, & Suhling, 2010; 
Norling, 1984b) and no accurate data on these thresholds are available for the studied species, a 
relatively broad period from 1 March to 30 September was included in the analyses. Two types of 
temperature variables were used. First, we calculated the mean water temperature for the period 
from 1 March to 30 September for both years separately (T_2: mean water temperature from 1 
March to 30 September two years before emergence; T_;: mean water temperature from 1 March to 
30 September one year before emergence). In addition, we used the temperature difference (AT_2, 
AT) between the actual mean temperature (T_2, T_;) and the long-term mean temperature (T20, 
encompassing 1989-2009) of the same period, i.e. 1 March to 30 September. 


Statistical analyses 


To reveal the relationship between the sex ratio and water temperature, percentages of males 
(%S, SEM 25, CEM 0-35, %CEM 75-109 and %7 EMoo-100) were used as dependent and 
temperature data (T_2, T_;, AT_2, AT_}) as independent variables in correlation tests. We used 
Pearson correlation, since all datasets had normal distribution (Jarque-Bera test, p > 0.05 in all 
cases). 

First, tests were run including data from Tiszafüred (n = 4) only, which are likely to represent 
more reliable sex ratio records because of the higher exuviae numbers, exceeding the recom- 
mended minimum of 99 exuviae. Thereafter, we ran the same test including all sample sites 
(n = 7). Although in this case some sex ratio records were based on fewer than 99 exuviae, 
these results still may contribute to revealing the main tendencies. In the latter tests we utilized 
only AT_z and AT_, temperature variables to exclude local temperature differences between the 
sampling sites. For statistical analyses the program package PAST (Hammer, Harper, & Ryan, 
2001) was used. 
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Table 3. Temperature characteristics of the study sites and years. 





Collection site Collection season T20 Tı AT: T_» AT_» 
Tiszafiired 2001 16.9 17.7 0.8 17.1 0.2 
2002 17.0 0.1 17.7 0.8 
2003 18.3 1.4 17.0 0.1 
2009 16.8 —0.1 18.5 1.6 
Jánd 2008 14.3 15.6 1.3 14.0 —0.3 
Vásárosnamény 2008 15.0 16.9 1.9 15.8 0.8 


Olcsva 2008 15.9 17.7 1.8 15.5 —0.4 


Notes: T29: mean water temperature from 1 March to 30 September based on 20-year dataset; T_;: mean water temperature from 
1 March to 30 September one year before emergence; AT_,: difference between T_; and T29; T_2: mean water temperature 
from 1 March to 30 September two years before emergence; AT_,: difference between T_2 and T29. 


Table 4. Results of Pearson correlation. R-values are recorded where correlation is significant (except for 
two cases, where strong although not significant correlations were found). 














Gomphus flavipes Gomphus vulgatissimus 
Pairs of variables Tiszafiired All sites Tiszafiired All sites 
%S -Tı 0.9574* — 0.94607 — 
%S —- AT_1 0.9574* 0.9086** 0.94647 —0.9268** 
%SEM75-100 — AT_-1 ns (0.7162) 0.8325* ns ns 
%SEM 0-100 — AT_1 ns (0.8198) 0.8161* ns ns 


Notes: %c": percentage of males of the total; %@EM75_100: percentage of males of the total between EM75 and EM100; 
%S%EMop-100: percentage of males of the total between EMoo and EMjo0; T_1: mean water temperature from 1 March to 
30 September one year before emergence; AT_,: difference between T_; and T29. 

*p < 0.05, **p < 0.01, +: marginally significant (0.05 < p < 0.06); ns: not significant, —: correlation was not analyzed 
(see “Materials and methods”). 


Results 


No obvious bias towards one sex could be detected in either species, but sex ratio varied between 
sites and years from a slight excess of females to slight excess of males in G. flavipes (48.0- 
53.2% males), and ranged from a somewhat higher excess of females to slight excess of males in 
G. vulgatissimus (43.2-52.0% males) (Table 2). 

T_, was slightly higher than T29 at all sites and in every year, and the difference was negative 
in only one case. Similarly, T_. was lower than T29 only in two cases (Table 3). In addition, no 
obvious trends could be observed in changes of AT_; and AT_2. 

T_2 and AT_» temperature variables were not correlated with %c’, but T_; and AT_; showed 
a significant relationship with %o’. These correlations were positive in G. flavipes but negative in 
G. vulgatissimus. Moreover, %8 EM75-100 and %o" EMoo_100 Were in strong positive correlation 
with AT_, in G. flavipes (Table 4), although the relationships were only significant when all sites 
were included in the analysis. 


Discussion 


In our study, neither G. flavipes nor G. vulgatissimus exhibited a clear bias towards one sex, but both 
had variable sex ratios. This result is consistent with some long-term studies that recorded no con- 
stant bias, but found great inter-year variability. Kern (1999) reported male ratios of 42.5-52.9% 
for G. vulgatissimus from a small stream, while Müller (1995) demonstrated male ratios of 30.7— 
72.4% for G. vulgatissimus and 40.1—51.0% for G. flavipes from the River Oder. In two short-term 
studies much more similar sex ratios were also recorded for G. vulgatissimus: Beutler (1986) 
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reported male ratios of 41.8% and 46.1% from the River Spree, while Westermann, Westermann, 
Heitz, and Heitz (1995) found male ratios of 50.0% and 50.2% at two different sites along the 
Upper Rhine. 

Despite the small dataset used in this study our results suggest that a relationship exists between 
water temperature in the year preceding emergence and the sex ratio at emergence. However, the 
way in which temperature affects the sex ratio is still unclear. Under natural conditions it is unlikely 
that water temperature can change to an extent that would be lethal to the larvae. Therefore, water 
temperature during larval development may have indirect effects on the sex ratio at emergence. 
Sex-related response to varying water temperature was described by De Block and Stoks (2003) 
but such responses were not examined in our study. 

However, a sex-specific effect of temperature is expected to be similar in the two closely 
related Gomphus species (with overlapping microdistribution, similar behavior and prey organ- 
isms; see Miiller, 1995). Contrary to this expectation the correlation between water temperature 
and sex ratios was of opposite direction in the two species, being negative in G. vulgatissimus, 
but positive in G. flavipes. Therefore, it is reasonable to suggest that the correlation might be 
attributable to other developmental characteristics. One possible explanation is related to voltin- 
ism and cohort-splitting, along with the faster development of males. Although this explanation 
is largely theoretical, it is supported by our results if the following elements are taken into 
consideration. 

It has been reported that cohort-splitting (a part of a year-class emerges either one year before 
[Corbet, 1957] or after [Norling, 1984a] the majority of that year-class) can influence the sex ratio; 
in Anax imperator a variably developed univoltine cohort late in the emergence period with female 
predominance affected the sex ratio at emergence (Corbet, 1957). Cohort-splitting is well docu- 
mented in gomphids (e.g. Aoki, 1999; Ferreras-Romero, Atienzar, & Corbet, 1999) and also in the 
studied species (Miiller, 1995). Among its causes, water temperature and intraspecific competition 
for food have been stressed (e.g. Banks & Thompson, 1987; Corbet, 1957; Norling, 1984a). Prey 
organisms, especially chironomids, are abundant in the River Tisza (e.g. Móra, Csépes, Toth, & 
Dévai, 2008), providing a nearly unlimited food resource for the gomphid larvae. Therefore, water 
temperature is more likely to be playing a major role in splitting year-classes in our study and the 
size of the minor cohort may vary year to year depending on the actual temperature regime. 

The unequal proportions of the sexes in the minor cohort may derive from the faster development 
of males than females: males that emerge earlier may have higher mating success, whereas it may 
be advantageous for females to develop longer, emerge at larger size and thus be more fecund 
(e.g. Corbet, 1999). If so, cohort-splitting may cause a 1-year earlier emerging cohort to contain a 
higher proportion of males, while a 1-year later emerging cohort contains a higher proportion of 
females. The faster development of male larvae has been demonstrated for a number of species by 
laboratory studies (Baker et al., 1992; Fuselier et al., 2007; Mikolajewski, Brodin, Johansson, & 
Joop, 2005), and may also be indicated by protandry, which is common in gomphids (e.g. Aoki, 
1999; Lutz & McMahan, 1973; Miiller, 1995; Suhling & Miiller, 1996; Westermann et al., 1995) 
and was exhibited in the populations studied here (see Jakab, 2006). 

The voltinism of Odonata is relatively flexible, depending on photoperiod and water temperature 
in temperate regions (Braune et al., 2008; Corbet, Suhling, & Soendgerath, 2006). At our sampling 
sites life cycles 2 to 3 years long are suggested for both studied species (Miiller, 1995; Miiller 
et al., 2000; Söndgerath et al., 2012). However, there is evidence that the development of G. 
vulgatissimus can vary greatly (from 1 to 4 years) between different habitats (Müller, 2000; 
S6ndgerath et al., 2012), while our knowledge on the life cycle of G. flavipes is more limited 
(Miiller, 1995). 

As has been demonstrated, larval development of G. vulgatissimus could be finished within two 
years in warm and shallow rivers (Müller et al., 2000; Sondgerath et al., 2012). Similarly, a 2-year 
life cycle was predicted for this species up to 50°N, based on modeling population dynamic along 
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Figure 1. Predicted life cycles of (a) G. flavipes and (b) G. vulgatissimus to explain the relationship between water 
temperature and sex ratio through cohort-splitting (t-2—t1: the years of development). 


a latitudinal gradient (Braune et al., 2008). In our study the rivers Tisza and Szamos met these 
criteria. In G. vulgatissimus, assuming a 2-year major and a 3-year minor cohort (Figure 1), the 
statistical correlations between sex ratios at emergence and the water temperature in the previous 
year (T_;, AT_;) may be explained as follows. After a year with higher water temperature (t—1), 
faster larval development can result in fewer individuals of G. vulgatissimus emerging one year 
later (tı) than the majority of that year class. This means that the 3-year minor cohort is relatively 
small, thus more females of the year-class emerge in their second year (to). This causes a decreased 
proportion of emerging males in tọ resulting in a negative statistical correlation between sex ratio 
and water temperature. After a year with lower water temperature (t_,), slower larval development 
can result in more individuals of G. vulgatissimus emerging one year later (tı) than the majority 
of that year class. Due to their faster development, relatively fewer G. vulgatissimus males than 
females of the year-class extend their larval stage to the third year (t,). This results a higher 
proportion of emerging males in to, thus supporting the negative statistical correlation between 
sex ratios and water temperature. Although it was not confirmed by the correlation tests, this 
mechanism can also lead to a negative correlation between water temperature and male ratio at 
the end of the emergence (%8 EM35-100 and %co" EMop_199). One reason that these correlations 
could not be detected may be the small numbers of exuviae. Moreover, for G. vulgatissimus in 
the River Oder, Miiller (1995) suggested a more complex scenario of cohort-splitting, with more 
than two cohorts emerging in the same season. 

In contrast, with G. flavipes having a similar life cycle the positive correlation between sex ratio 
and water temperature cannot be explained in the previous way, suggesting that the development of 
this species is different from that of G. vulgatissimus. Assuming a 3-year major and a 2-year minor 
cohort (Figure 1) however, the results can be easily interpreted. After a year with higher water 
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temperature (t_,), faster larval development can result in more individuals of G. flavipes emerging 
one year earlier (to) than the majority of that year class. Accordingly, the 2-year minor cohort is 
more developed. In addition, males develop faster than females, meaning that relatively more males 
of the year-class emerge in their second year (tg) and relatively more females of that year-class 
emerge in their third year (tı). Due to this an increased proportion of emerging males appears in to, 
resulting in a positive statistical correlation between the sex ratio and water temperature. After a 
year with lower water temperature (t_,), slower larval development can result in fewer individuals 
of G. flavipes emerging one year earlier (tg) than the majority of that year class. Thus, the 2-year 
minor cohort is less developed, which means that more males of the year-class emerge in their 
third year (tı). This decreases the proportion of emerging males in to, contributing to the positive 
statistical correlation between sex ratio and water temperature. Moreover, the 2-year minor cohort 
is liable to emerge later in season than the majority, which was confirmed by the correlation tests 
(positive correlation between water temperature and %0 EM75_109 and Yoo" EMoo-100). 

Additionally, in this explanation t_2 corresponds to the year of hatching in both species 
(Figure 1). In the year of hatching larval growth is restricted to a relatively short period, and 
therefore the water temperature of this year may have less influence on the development com- 
pared to that of the following year. This could explain why T_, and AT_, (the water temperature 
in the year of hatching) were not correlated with sex ratio in the statistical tests. 

Summarizing the relationships discussed above, the effect of water temperature on the sex 
ratio at emergence may be mediated through cohort-splitting. A variably developed minor cohort, 
including unequal proportions of males and females, may affect the sex ratio, where the size of 
this minor cohort depends on the water temperature. Note that a higher proportion of males does 
not necessarily exceed 50% in these cases, it simply means an increase in proportion compared 
to the preceding or subsequent year. Hence we can conclude that the effect of water temperature 
may not cause a long-term consistent bias, but might explain the year-to-year variations in sex 
ratio. 

The various studies on sex ratio of Odonata at emergence imply that mechanisms acting on sex 
ratio are complex (see “Introduction’”’). Nevertheless, in accordance with our findings, in future 
studies it would be reasonable to take into account water temperature, voltinism and cohort- 
splitting, as potential causes of the variability in sex ratios. 
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